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ABSTRACT 

For a long time globular clusters (GCs) have been considered as ideal laboratories to study dense stellar 
systems. However, recent observations have revealed that most of them consist of different stellar populations. 
In particular, omega Centauri (ui Cen) is one of the most complex GCs and has attracted a lot of attention, 
because its high mass may depict a link between GCs and larger stellar systems, such as ultra compact dwarf 
galaxies (UCDs). Here we show with direct-summation A/-body simulations that a natural explanation is the 
merger of GCs in colliding galaxies, such as the Antennae galaxy. Hubble Space Telescope observations of 
these interacting galaxies show clusters of clusters, regions in the galaxy that span hundreds of parsecs, where 
many of the GCs are doomed to collide. Also, observations of the Large Magellanic Clouds (LMC) reveal that 
a fraction of GCs are bound and thus will result in a merger. Our results prove that colliding GCs with different 
metallicities and ages can produce a GC with multiplicity and occupation fractions similar as observed. In our 
scenario the merged clusters have a phase with a larger amount of flattening than average, and we find that 
observational data correlates evidence of rotation with multiplicity, which corroborates our hypothesis. Future 
observations of multiple populations should focus on GCs with ellipticity above average. 



1. MOTIVATION 

Thanks to the Hubble Space Telescope (HST) we have a 
growing number of observations that contradict the traditional 
narrative about GCs being perfect laboratories for studying 
the structure and evolution of stars; i.e. simple stellar pop- 
ulations that are coeval and chemically homogeneous. Inde- 
pendent observations corroborate that a number of GCs har- 
bor stars which have formed from material from different 
generation of stars (see for a review Piotto 2009 and refer- 
ences therein). The first discove ry of multiple main sequences 
(MS), u Cen ( |Bedin et al.|2004) , turned out to be one of the 
most com plex ones, c ontaining not only two but a third, red- 
der MS ( |Villanova et al.|[2 007) with the blue MS with twice 
the metal abund ance of the MS dominant red branch ( |Pi-| 
lotto et al.|[2005l >. Other cas es followed; the CMP of NGC 
2808 showed a trimodal MS ( |Carretta et al.|2006[ [Piotto et al.| 
2007), NG C 1851 possesses at least two dissimilar stellar pop- 
ulations ( |Milone et al.|2008"| which manifest at the splitting of 
the sub-giant branch (SGB) in the CMD. In the case of Galac- 
tic GCs, we know of many other clusters with a bimodal SGB, 
as e.g . M22 ( jMilone et al.|20TT) , NGC 6388 and M54 (Piotto 
2009). In the case of the Large Mage llanic Cloud (LMC) up 
to 70% of the GCs show multiplicity ( |Milone et al.|2009a) . 

On the other hand, the merger of stellar clusters is a plau- 
sible situation in nature. Observations of colliding galaxies 
such as the Antennae galaxy show bound systems of young, 
massive clusters. In this system HST observations exhibit 
relatively small regions spanning a few hundreds of pc em- 
bracing hundreds or even thousands of young clusters, i.e. 
clusters of clust ers or "clu ster com plexes" (CC fro m now on- 
wards, see e.g. |Kroupa| 19 98; Whitmor e et al.|2010| and ref- 
erences therein). These have been proposed to be the progen- 
itors of UCDs as a result of the agglomeration of hundreds of 
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their member clusters (|Fellhauer & K roupa 2002 , 2005 ; Briins| 

1 lb . Also, in the Magel- 
lanic Clouds (MCs) a large fraction of clusters are observed 



et al.|201l"t |Amaro-Seoane et al.||2011|i. Also, in the Mage 



Dieb aU et al.| f2002) 



in bound binaries or low -order multiples, 
estimate that about 1 /8 of the clusters in the Large Magellanic 
Cloud is a member of a bound group. 

In this Letter we show that CCs and the LMC are a natural 
breeding ground for multiple-populations GCs because they 
collide and merge. We run direct-summation A/-body simu- 
lations and analyse the impact of dynamics on the occupation 
fractions of the different populations in the CMDs. We show 
that two colliding GCs produce a system with phases of el- 
lipticity significantly over average and, as a matter of fact, 
observational data reveal that GCs harboring multiple stellar 
populations are likely to have a larger amount of rotation or 
ellipticity than average values, so that future surveys of mul- 
tiple stellar populations should focus on these. 

2. SMASHING CLUSTERS 

We set initially the clusters on a parabolic orbit so that the 
minimum distance at which they pass by is of m j n if they are 
considered to be po int particles at their centers of mass, as 
described in |Amaro-Seoane ( 2006| l. We integrate the evolu- 
tion with direct-summation A^-body tools, which integrate all 
gravitational forces for all particles at every time step, without 
making any a priori assumptions about the system. The codes 
we have employed, A/-body a nd Myriad use the improved 
Hermite integration sche me ( |Aarseth||1999| |2003| [Konstan- 
tinidis & Kokkotas 2010). This needs computation of not only 
the accelerations, but also their time derivatives. The pro- 
grammes also include KS regularisation and chain regulari- 
sation, so that when particles are tightly bound or their sepa- 
ration becomes too small duri ng a hyperbolic encounter, the 
system is regularised (Kustaanheimo & Stief el| 1965} |Aarseth] 
2003) to prevent too small individual time steps. Direct- 



summation A/-body codes scale as N+, or Af cx td yn , with fdyn 
the dynamical time. For this reason we have employed a ver- 
sion of Myriad which runs on the GRAP E, and the GPU 
version of A/-body of Sverre Aarseth ( Aarseth 2003 ; Ko nstan-| 



2 





Parameters 






ap"P ( M vv) 
\ lvi y L f 


Distances 


Case A 


Cluster 1 


9 


0.002 


100 


^COM = 12 
rfmin = 0.5 


Cluster 2 


5 


0.001 


50 


CaseB 


Cluster 1 


12 


0.006 


50 


«COM = 10.6 

*^min — 0-5 


Cluster 2 


5 


0.005 


100 


CaseC 


Cluster 1 


12 


0.02 


50 


AtoM = 10.6 

dmm = 0.5 


Cluster 2 




0.01 


100 


Case© 


Cluster 1 


9 


0.002 


100 


^COM = 52.8 
dais = 20 


Cluster 2 


5 


0.001 


10 


Case £ 


Cluster 1 


9 


0.002 


100 


1+2 from .A 
^COM,2=49 
4nin.2 = 10 


Cluster 2 




0.001 


50 


Cluster 3 


5 


0.009 


2.1 x 10* 



TABLE 1 

Initial conditions for the clusters. From the left to the 
right, we have the king parameter, initial metallicity and 
the initial age of the cluster in myrs. on the last right 
column we show the initial distance between the 

CENTER-OF-MASS OF THE TWO CLUSTERS THA T COLLIDE, Rroivi 
(WH ICH CORRESPONDS TO |xi+X 2 | OFFlG. 1 IN |AMAR0-SE0ANE| 
( |2006| >), AND d min , BOTH IN PC. IN CASE £ WE USE THE OUTCOME OF A 
AND MAKE IT COLLIDE WITH A THIRD CLUSTER, OF N* = 20,000 STARS 
FOR £ . IN THIS SIMULATION WE ADD THE NEW DISTANCE BETWEEN 
COM, Rqom.2 AND fif min .2, FOR THE SECOND COLLISION. 

|tinidis & Kokkotas|2010) . 

3. CMDS: MULTIPLE POPULATIONS AND FRACTIONAL 
OCCUPATION NUMBERS 

We ran different cases with different initial conditions, as 
shown in table [1] The clusters were modelled initially with a 
King model of different concentrations and metallicities and 
evolved for differ ent times with the stellar evolution package 

(2000). For the particular case 



sse, described in Hurley et al. 



£ , we use the outcome of simulation A and make it merge 
with another cluster (numbered 3) on a parabolic orbit with a 
new c/ m i n .2 and /?com,2> as indicated on the right of the table [T] 

3 During the collision we neglect stellar evolution because in 
1 runs it took approximately a few Myrs and the impact of 
evolving the masses on the global dynamics is negligible. We 
find in our simulations a significant mass loss after the merger 
of the clusters that affects the different occupational fraction 
numbers of the CMDs. 

To determine when the clusters me rge, we locate the density 
centers of the two clusters following Casertano & Hut ( 1985 ) 
and that of the merged system. We stop the simulations when 
these three centers have coincided for a few fd yn . This allows 
us to study the distribution of stars due to the dynamics of 
the system, which is important to understand the impact of 
mass loss in the CMD and occupational fractions in different 
shells of mass around the density center of the merged sys- 
tem. The CMD we obtain from the simulations correspond 
to idealised observational conditions. In real observations the 
measurements are affected by systematic errors. For instance, 
in the figure 10 of Bellini et al.| ( |2010| , which corresponds to 
the calibration field of w Cen fortEeHST WFC3/UVIS with 
the filters F275W, F814W, the first one has a RMS of sC 0.01 
in magnitude for well-measured, unsaturated stars. Saturation 
for F814W starts around the horizontal branch (HB) level and 
can increase to ~ 0.05-0.08 magnitude for faint stars in the 
CMD, while saturated stars in F814W can have errors up to 
0.05 in magnitud^] We use this as a representative exam- 
ple and introduce an artificial error in our CMD in F606W 



4 For movies see, |http : //www, aei .mpg. de/~pau/| 
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of ±0.01, and ±0.05 in F814W by scattering the data with 
a Gaussian distribution around the real data to show how the 
populations would appear in an observation of the system. In 
Fig.([T]i we show the CMD of simulation £ . We can see that at 
the level of the turn-off point (TO) the SGB splits into three 
branches. 

We then calculate the occupation ratio of the different stel- 
lar populations, N\ /N2 (and Ni/[N\ +N2] if we had three dif- 
ferent populations) for different shells starting at the density 
center. The results are shown in table [2] From an obser- 
vational point of view, there is a debate on whether or not 
there is a significant difference in the distribution of the dif- 
ferent populatio ns in GCs wit h multiplicity. In the case of 
NGC 1851, |Carretta et al. | pOlO) [2011) observed evidence of 
a difference in the distribution of the two populations of red 
giants. According to them, the metal-poor (MP) component 
seems to be more centrally concentrated than the metal rich 
(MR) one. This observational evidence contradicts the theory 
of the formatio n of the se cond population of stars within the 
same cluster (Bekki 201 1), in which the metal rich population 
should be more concentrated around the center. On the other 
hand, [Milone et alT| | |2009b) did not find any significant radial 
trend in the number ratio b etween the two p opulations of stars 
of the SGB (although see |Carretta et al.| < 



2011) . The distri- 
bution of the different populations in the final cluster of our 
models depends mostly on the initial concentration, metallic- 
ities, and also on the initial number of stars in each cluster, 
as well as their initial mass function and ages. In most of our 
models (all but B and C), the metal-rich cluster, which is al- 
ways cluster 1, appears to have fewer stars in the center of the 
final cluster than the metal-poor population, which is a natu- 
ral result of the initial King parameter. u> Cen is the the sys- 
tem for which we have the best observational data about the 
radial distribution of its multiple populations. Belli ni et al.| 
(2009) do a detailed study showing that stars belonging to the 
blue MS (bMS) appear to be more centrally concentrated than 
stars of the red MS (rMS), with the fraction A^ms/Mms < 1 
outside the core of the system. Since the bMS contains stars 
with greater metallicity, according to the authors, stars with 
greater metallicity appear to be more centrally concentrated. 
As shown in table [2j one of our models (B) shows a distri- 
bution in which the MR stars appear more centrally concen- 
trated relative to the MP stars of the final cluster. In the same 
model, metal poor stars appear to dominate the innermost re- 
gion, while metal rich stars again dominate the external shells 
of the cluster as a result of the large difference in the initial 
concentration of the two clusters. 

4. RESULTS AND CONCLUSIONS: THE ROLE OF ELLIPTICITY AS 
A SIGNATURE 

In this paper we have addressed the observational evidence 
of multiplicity of stellar populations in globular clusters in 
terms of collisions of clusters. Interacting galaxies such as the 
Antennae or the LMC are natural loci for clusters to collide. In 
the later, CCs have been observed with the HST and they are 
the natural birth -place of UCDs, as proposed by a number of 
works dKroupa|1998||Fellhauer & Krou pa|2002"l |20"05| [BTuns] 
|et al.|201 l[|Amaro-Seoane et al.|201 1) . 

While dynamics does not affect the shape of the resulting 
CMD of the merged system, it impinges the number of stars in 
different parts of it: the occupation fraction will vary because 
of the dynamics. 

After the collision, the resulting cluster has a significant 
amount of rotation. This depends on the initial conditions 
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Ay/V 2 (andAVM+Ay) 


Shell (pc) 


Case A 


Case li 


Case C 


Case IJ 


Case £ 


Osg r<0.5 


0.94 


1.25 


1.08 


0.50 


0.85 (0.32) 


0.5 < r< 1 


0.86 


1.02 


1.03 


0.53 


0.80 (0.26) 


1 ^ r<2 


0.72 


0.83 


0.87 


0.61 


0.78 (0.23) 


2^r<3 


0.66 


0.81 


0.82 


0.90 


0.74 (0.23) 


3 r<4 


0.82 


1.00 


0.98 


1.07 


0.79 (0.24) 


4sJ r<5 


1.03 


1.06 


1.00 


1.25 


0.77 (0.27) 


5^r< 10 


1.39 


1.19 


1.17 


1.52 


0.95 (0.30) 


10^r<50 


1.96 


1.42 


1.36 


2.28 


1.31 (0.45) 


e 5 . max 


0.178 


0.215 


0.100 


0.174 


0.298 



TABLE 2 

Occupation fraction and maximum and maximum ellipticity 

FOR THE MASS FRACTION 5, e5. raax , FOR THE CASES OF TABLeR] THE 
FRACTIONS ARE GIVEN IN TERMS OF NUMBERS OF STARS BELONGING 

INITIALLY TO THE FIRST CLUSTER (N\) AND THE SECOND ONE (N 2 ) FOR 
DIFFERENT SHELLS OF THE RESULTING MERGED CLUSTER STARTING 

FROM THE DENSITY CENTER. FOR THE LAST TWO CASES WE ALSO GIVE 
THE FRACTION RELATIVE TO THE THIRD ONE, N 3 . 




0.4 0.6 0.8 1.0 

m F605W — m F814W 

FIG . 1 . — CMD of case £, a collision of a cluster with another cluster which 
is the result of a merger of two clusters itself (case A) with artificial errors 
(upper panel) and without them (lower panel). The first rMS (curve with the 
TO point at~ 0.45 in OJF606W — T814 m the lower panel, blue line in the on- 
line version of the article) corresponds to the stellar population of z = 0.002, 
the next rMS to z = 0.001 (red curve of lower panel) and the third one to 
Z = 0.009 (green curve). The clusters that harboured the two first metallicities 
had an age of 100 and 50 Myrs respectively, and the third one had an age of 
10 3 Myrs. The CMD is for stars in the shell ^ r/pc < 2.5 after evolving 
the merged cluster for 9 Gyrs. We set the distance of 4850 pc to convert to 
apparent magnitudes. 



such as the impact parameter, the King parameter and the rel- 
ative velocity. A comprehensive study of the parameter space 
and impact on the oblateness of the merged systems as well as 
the evolution of rotation and CMDs is out of the goals of this 
work and will be presented elsewhere. Indeed, we have di- 
rect measurement of the rotational velocity of a few GCs that 
corroborates a correlation between rotation and multiplicity 
of stellar popu lations. In most of the o bservations, rotation 
is differential (Meylan & Mayor 1986). The results for the 
rotati onal veloci ties of 11 GCs are summarized in Table 7.2 
of |Meylan & Heggie| ( |1997 1, where the rotational velocities 
v rot are measured in terms or the velocity dispersion a of each 
cluster. In Table [3] we show the e and the multiplicity of a 
few GCs for which we have data: 9 out of the 11 clusters 
show evidence of multiplicity. For the other two there are no 
observational data about multiplicity to our knowledge. The 



Cluster 



NGC5139 
NGC 6656 
NGC 3201 
NGC 7089 
NGC 6205 
NGC 6341 
NGC 104 
NGC 6397 
NGC 7078 
NGC 5272 
NGC 362 



0.17 
0.14 
0.12 
0.11 
0.11 
0.10 
0.09 
0.07 
0.05 
0.04 
0.01 



Multiplicity 



0.32 
0.50 
0.28 
0.34 
0.25 
0.30 
0.26 
0.11 
0.15 
0.12 
0.01 



3MS,6SGB,4RGB (Belli ni et al.|2010} 
2SGB.2R GB jM ilone et al. 20TT> 
2SGB,2RGB |Kravtsov et al.|2 010 
N o data available 
3RGB I Carretta et al. 



2011 



2009 



No data available 



2SGB,2RGB,broad MS I Anderson et al. 



2RGB jCarretta et al.pdof} " 
2RGB (Carretta et al.|2009| 
2RGB (Carretta et al.|2009| 
2RGB (Pancinoetal. 20101 



2009 



TABLE 3 

ELLIPTICITY, ROTATIONAL VELOCITY AND MULTIPLICITY WITH 
REFERENCE OF 11 GALACTIC GCS. 




60 65 

T (Myrs) 

FIG. 2. — Evolution of e for case T> after the density centers coincide. 
The semi-major axes are calculated with the ellipsoids of inertia (see e.g. 
Chandrasekhar 1969 1 and determined by different mass fractions of the stars 
(at 1 ~ 47.7 Myrs, from the bottom to the top 0.2, 0.5 and 0.8, displayed 
in orange, red and blue, respectively in the on-line article). The solid black 
line shows e; thus, values within the colored box are average. The stars are 
distributed according to the amount of gravitational energy; hence, the lower 
the mass fraction is, the closer we are to the centre of the resulted merged 
system. 

fastest rotating cluster is NGC 6656 with v m t/u = 0.5, and an 
ellipticity of e = 0.14 (with e := l-b/a, where a and b are 
the semi-major and -minor axes), to Cen rotates very fast, 
Vrot/°' = 0.32-0.41, which makes it one of the fastest rotat- 
ing GCs in the Galaxy. This could be a signature for an ag- 
glomeration process of a cluster in a CC which receives more 
and more impacts from other lighter clu sters and "runs away" 
in ma ss, on its way to forming a UCD ( |Amaro-Seoane et al.| 
[20TTT >. 

In Fig. pi we show e vs time for T>. We start the analysis 
after the tnree density centers coincide. The system has an e 
above average during some few Myrs and by T ~ 75 Myrs 
e ~ .08. We have data for e for 100 out of the 157 Galactic 
GCs ( |Harris|[2qT0l >; the average is e - 0.0766. Some of the 
clusters for which there is observational evidence of the ex- 
istence of multiple stellar populations have larger values than 
e. In particular, lu Cen has e = 0.17, which makes this clus- 
ter 9th in the list of the most oblate clusters. Other clusters 
have a rather large value, NGC 2808 (e = 0.12), NGC 6656 
(e = 0.14) and NGC 3201 (e = 0.12). 47 Tuc, with multiple 
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stellar populations (Anderson et al. 2009 1, has a value slightly 
above average, 0.09. 

We also note that there is no observable difference in the 
rotation of the multiple populations in G Cs. Measurements of 
the rotation of the 650 stars of uj Cen ( |Pancino et al.||2007) 
Anderson & van der Marel 2010) show that all subpopula- 
tions rotate as a single one. This can be explained in terms 
of a collision between GCs, which would assign all stars the 
same amoun t of rotation regard less of their population. Also, 
according t o|Lane et al.| p010) several of the observed prop- 
erties of 47 Tuc can be explained if a merger of two clusters 
occurred in the past. |Kravtsov et al.| ( |201 1 J > find that there is 
a strong segregation between the two populations of stars that 
exist in NGC 6752 and NGC 3201, with the fainter SGB and 
redder RGB stars dominating the center of both clusters. This 
difference in the distribution of the two populations might be 
the result of a merger between two GCs, as shown by our sim- 
ulations. 

Hence, although this mechanism probably cannot explaint 
all cases observed, in particular the strange case of the metal- 
licities of the different MS of uj Cen, we deem it is a plausible 
and natural way of creating multiplicity of stellar populations. 
Also, even if observations on rotation are incomplete and do 
not cover all of the GCs with multiplicity, we have shown in 
our limited study of the parameter space that the result of a 
collision is a cluster that will exhibit phases in the evolution 
with e above average. In the particular case of simulation T), 
the GC achieves average values after a short time, so that any 
oblateness would not be present today in old clusters, unless 



the collision happened recently (though we note that cannot 
model realistic GCs with our number of stars). This means 
that GCs, in particular young ones, with e above average are 
more likely to harbor multiple populations of stars; i.e. any 
amount of rotation in GCs with multiplicity could be a fin- 
gerprint for a dynamical origin. Future observations should 
focus on such GCs, since these are likely to display multiplic- 
ity. In particular, the systems NGC 7089 and NGC 6341 have 
not been investigated for presence of multiplicity, and in our 
scenario they are good candidates for having different stellar 
populations. 
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